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I. INTRODUCTION 

Motivated from string theory, D-brane physics and 
black hole thermodynamics, we focus on investigating 
the possibilities of covariant thermodynamic geometric 
study of hot QCD near the QCD transition tempera- 
ture. In particular, we intend to study the thermody- 
namic geometry arising from the free energy for the 2- 
and 3-fiavor finite temperature QCD. The main idea is 
first to develop a geometric notion for QCD thermody- 
namics and then to relate it with the existing microscopic 
quantities already known from hot QCD [J 0. We fur- 
ther incorporate thermal fluctuations in the free energy 
[H and thus study their contributions to the covariant 
thermodynamic geometry. 

Let us recall that the existence of a geometric struc- 
ture in equilibrium thermodynamics was introduced by 
Weinhold [4| through an inner product in the space of 
equilibrium thermodynamic macrostates defined by the 
minima of the internal energy function U = U(ui,v,s) 
as the Hessian function — didjU. Here the reduced 
quantities {vi,v, s} may respectively be defined as un- 
derlying chemical potentials fa, volume V and that of 
the entropy S scaled with temperature of the considered 
equilibrium configuration. In order to provide a physical 
scale and to restrict negative eigenvalues of the metric 
tensor, we require the volume or any other parameters to 
be held fixed for a given QCD configuration. 

Interestingly, the Riemannian geometric structure thus 
involved indicates certain physical significance ascribed 
to hot QCD. The associated inner product structure on 
the intrinsic space may either be formulated in the en- 
tropy representation, as a negative Hessian matrix of the 
entropy with respect to the extensive variables, or that of 
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the energy representation, as the Hessian matrix of the 
associated energy with respect to the intensive variables. 
Such geometric considerations have been studied earlier 
in the literature as the covariant thermodynamics whose 
application in the energy representation has first been 
considered by Weinhold 4]. On the other hand, the en- 
tropy representation has been considered in Ruppeiner's 
fluctuation theory Q . 

In this paper, we shall take the former picture and 
study the covariant thermodynamic geometric properties 
ascribed to hot QCD. The thermodynamic macrostates 
of an underlying equilibrium chemical potential space 
may in turn be described by the minima of the energy 
function F = F(U,V,N). Explicitly, the thermody- 
namic metric tensor of an intrinsic Riemannian space 
spanned by the chemical potentials may be given as 
9ij = 9idjF(fj,i,T, V, S) and turn out to be conformal to 
the Ruppenier metric with the temperature as the con- 
formal factor. This intrinsic metric as an inner product 
on the (hyper-) surface of chemical potentials is thus nec- 
essarily ensured to be symmetric, positive definite and 
satisfies the triangle inequality since the energy function 
has a minimum configuration in the equilibrium. 

On other hand, the thermodynamic geometry of the 
equilibrium configurations thus described has extensively 
been applied to study the thermodynamics of a class of 
black hole configurations [7|. Recent studies of the ther- 
modynamics of diverse black holes in the state-space ge- 
ometric framework have elucidated interesting aspects of 
phase transitions and their relations associated with the 
extremal black hole solutions in the context of the mod- 
uli spaces of N > 2 supergravity compactifications [1]. 
It may be argued however that the connection of such 
geometric formulations to the fluctuation theory of black 
holes thermodynamics requires several modifications 0] . 
The geometric formulation involved has first been applied 
to N > 2 supergravity extremal black holes in D = 4 
which arise as low energy effective field theories from 
compactification of Type II string theories on Calabi- 
Yau manifolds [Tc|. Since then, several authors have 
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attempted to understand this connection [llj, Il2| both 
for supersymmetric as well as non-supersymmetric black 
holes and five dimensional rotating black rings. Such geo- 
metric notions further remain well-defined for the diverse 
cases of extremal as well as non-extremal black branes, 
see for details @- 

Here, we shall study the standard two dimensional 
thermodynamic geometry associated with the space of 
chemical potentials {fJ>i, ^2} in the framework of quasi- 
particle theories of QCD. The application of our geo- 
metric notions, as in any conventional thcrmodynamical 
systems, suggests that the space defined by the chemical 
potentials indicates the nature of the underlying statisti- 
cal system. In particular, such geometric notions like the 
covariant nature of the metric tensor and non-zero scalar 
curvature may shed light on the interactions present in 
the system. It may further be shown that the confor- 
mally related thermodynamic scalar curvature actually 
defines the associated correlation length of any arbitrary 
dimensional physical system. 

As a matter of fact, a few simple manipulations il- 
lustrate that the conformal thermodynamic curvature is 
inversely proportional to the singular part of the free en- 
ergy associated with long range correlation(s) whose di- 
vergences consequently signify certain phase transition(s) 
at the critical point(s) of the QCD. It is worth to mention 
that the Ruppeiner formalism may further be applied to 
diverse condensed matter systems with two dimensional 
intrinsic Ricmannian spaces, which has been found to be 
completely consistent with the scaling and hyper scal- 
ing relations involving certain critical phenomena and in 
turn reproduces the corresponding critical indices, see for 
details [5[. 

The purpose of this paper is to investigate the rela- 
tion between the covariant thermodynamic geometry and 
the quark-number susceptibility tensor, both of which 
have recently received considerable attention from the 
perspective of black hole physics and QCD thermody- 
namics. Thus the goal here is to explore underlying mi- 
croscopic arguments behind the notion of the covariant 
thermodynamic geometry arising from the free energy 
of certain QCD thermodynamic configurations. It has 
however been known in quasi-particle models that the 
quark number susceptibilities may in general be defined 
7tjr~ = fh,. an. ~ Xjt- Here i, j are flavor indices, 



Xij 



Mi is the quark number density, and P is the associated 
pressure of the QCD configurations, see [l[ for details. 

In this framework, it has further been shown that 
lattice measurements of the off-diagonal susceptibility 
Xud determine the associated physical properties of the 
quark-gluon plasma. As reported in Ref. Q, these are 
shown to be consistent with Hard Thermal Loop (HTL) 
studies down to the temperature of order as low as that 
of 2.5Tc, see also 13|. Blaizot, Iancu and Rebhan have 



in particular analytically computed the quark number 
susceptibilities for the quark gluon plasma at finite tem- 
perature within the self-consistent resummation of Hard 
Thermal/Dense Loop (HTL/HDL) [14]. Interestingly, 



these predictions at zero chemical potential were found 
to be consistent with the lattice results. 

Note that all such diagonal and off-diagonal elements 
of the Xij become equal at fii = 0, and thus the quark 
number susceptibility tensor reduces to two different sec- 



tors, and in particular, we have Xij 



and that of Xij 
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X for 



J, 



= x f° r i j- We have shown 

that the same remains true for the components of the co- 
variant thermodynamic metric tensor which may further 
be regarded as the above quark number susceptibility 
tensor. This equivalence in turn explains the involved 
macroscopic/ microscopic duality relations in our geo- 
metric investigations of the 2 and 3-flavor QCD near tran- 
sition temperature and turns out to be in close connec- 
tion with the AdS/CFT correspondence [HI]. Our anal- 
ysis thus clarifies the duality relationship between the 
macroscopic thermodynamic geometry and the micro- 
scopic quark number susceptibility formalism away from 
the ^ ^ 0. We note further that the microscopic formal- 
ism is based on the evaluation of the diagonal susceptibil- 
ity x within a chosen resummation scheme which directly 
follows from the expression for the fermion number den- 
sity Af in terms of the associated dressed fermion prop- 
agators as a function of the self-energies and fermions 
helicities, see for details [ill ]. 

On other hand, off-diagonal susceptibility tensors may 
be classified by the symmetry under charge conjugation, 
or that of C-parity. In particular, this in terms of chem- 
ical potentials means that a quark loop with two gluon 
external lines is necessarily even in the chemical poten- 
tials, while a quark loop with three gluon legs may also 
generate a linear term in the chemical potentials, which 
turns out to be symmetric in the color indices. We thus 
find a perfect match that our even derivative geomet- 
ric quantities like the components of the metric tensor 
or the associated determinants and the thermodynamic 
scalar curvature as the fourth derivative invariant quan- 
tity turn out to be even in the chemical potentials. We 
envisage that the components of thermodynamic metric 
tensor microscopically correspond to those of the suscep- 
tibility tensor and thus require two fermion loops con- 
nected by certain gluon lines. 

As the final exercise, we shall discuss the thermody- 
namic geometry of hot QCD including thermal fluctu- 
ations about the equilibrium configuration of chemical 
potentials. As it is well known, any thermodynamical 
system which may be considered as an ensemble has log- 
arithmic and polynomial contributions to the free energy 
Q. Moreover, this connection have earlier been explored 
from the perspective of black hole physics, showing that 
the thermodynamic Ruppeiner geometry of non-extremal 
rotating BTZ black holes, as well as that of BTZ-Chern 
Simons black holes, get modified under such thermal fluc- 
tuations. See for example 17[ , in which it has been shown 
that the thermodynamic geometry is flat for both such 
black holes, with or without certain higher derivative con- 
tributions. 
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What follows here, within the context of QCD, is that 
such geometric notions implicitly assume a statistical ba- 
sis in terms of an ensemble, and thus may explore the 
role of thermal fluctuations in the free energy and conse- 
quently on an associated intrinsic Riemannian geometry 
corresponding to quasi-particle theories. For example, 
the canonical free energy corrected by the contributions 
coming from the thermal fluctuations may be regarded 
as a closer approximation to that of the free energy in 
the corresponding ensemble. We here show that such 
considerations apply as well to hot QCD considered as 
an ensemble, and the specific forms of the logarithmic 
and polynomial contributions may easily be calculated 
for a wide class of quasi-particles theories. It is to be 
noted that the applicability of this analysis presupposes 
that the underlying ensemble is thermodynamically sta- 
ble, which requires a positive specific heat, or correspond- 
ingly that the Hessian matrix of the free energy function 
must be positive definite. The free energy function for 
any thermodynamical system incorporating such contri- 
butions @ may easily be written to be F = F +ln(CT 2 ). 
In this article, we have shown that such logarithmic con- 
tributions to the free energy induces bumps in the scalar 
curvature of the intrinsic thermodynamic space spanned 
by the quark chemical potentials. 

The paper is organized as follows. In the first section, 
we have presented the motivations to study the intrinsic 
Riemannian space obtained from the Gaussian fluctua- 
tions of the QCD free energy. In particular, we have 
outlined some microscopic implications of the geometric 
QCD in the lieu of the AdS/CFT correspondence and the 
remaining significance has been summarized in the final 
section. In section 2, we have introduced very briefly 
the quasi-particle thermodynamic geometry. In section 
3, we obtain the 2- and 3-flavor QCD free energy and 
temperature near the Tq as a function of quark chemi- 
cal potentials and then incorporate thermal contributions 
into the free energy. In section 4, we analyze the thermo- 
dynamic geometry of two-flavor QCD in the framework 
of the quasi-particle theory with or without logarithmic 
contributions. 



II. THERMODYNAMIC GEOMETRY 

In this section, we present a brief review of the essential 
features of thermodynamic geometries from the perspec- 
tive of the application to the hot QCD. In particular, 
we shall focus our attention on the geometric nature of 
the QCD thermodynamics in the neighborhood of QCD 
transition temperature, Tc with a finite number of quark 
chemical potentials introduced within the framework of 
quasi-particles theories. 

Let us consider an intrinsic Riemannian geometric 
model whose covariant metric tensor may be defined as 
the Hessian matrix of the free energy with respect to 
an arbitrary finite number of chemical potentials carried 
by the quasi-particles in a given thermodynamic config- 
uration at a fixed volume or any other parameters of 
the equilibrium QCD. In particular, let us define a rep- 
resentation F(fii,T) for given any free energy, chemical 
potentials, and temperature {F, /ij,T} and then choose 
a particular temperature slice as an intersection of the 
line T := with the space spanned by the underlying 
chemical potentials {/i^, T}, that is to define the temper- 
ature as in Eq. [TTJ In this consideration, the free energy 
thus defined may further be shown to be only a function 
of the quark chemical potentials, which in fact we have in- 
dicated in general by Eqs. [T^IHH This method in general 
turns out to be an intrinsic Riemannian geometry, which 
is solely based on nothing other than the quark chemical 
potentials. This thus yields that the space spanned by 
the n chemical potentials of the theory under considera- 
tion exhibits a n-dimcnsional intrinsic Riemannian mani- 
fold M n . The components of the covariant metric tensor, 
i.e. the so called thermodynamic geometry [j| @, [12, G3 
may be defined as 



where the vector x — (//) £ M n . Explicitly for the 
case of two dimensional intrinsic geometry, the com- 
ponents of the thermodynamic metric tensor are given by 



In section 5, we focus our attention on the thermody- 
namic geometry of the 3-flavor QCD configurations. Ex- 
plicitly, we investigate the nature of the thermodynamic 
geometry thus defined as an intrinsic Riemannian mani- 
fold obtained from the free energies with or without the 
contributions being considered in the section 4. We have 
explained that such a geometry obtained from the free 
energy of the quasi-particle theories results to be well- 
defined and pertains to an interacting statistical system. 
Finally, section 6 contains certain concluding issues and a 
discussion of the thermodynamic geometry of both the 2- 
flavor and the 3-flavor hot QCD and further microscopic 
implications arising from the underlying QCD, and thus 
offers a physical explanation of finite chemical potential 
macroscopic/microscopic duality relations in close con- 
nection with the AdS/CFT correspondence. 
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Rij and R for above 



Now we can calculate the , Rijki 
two dimensional thermodynamic geometry (M2,g) and 
may easily see that the scalar curvature is given by 
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Furthermore, the relation between this thermodynamic 
scalar curvature and the thermodynamic curvature ten- 
sor for such a two dimensional thermodynamic intrinsic 
geometry is given (see for details) [jj by 



R -y\\ R 



(4) 



which is quite usual for any intrinsic Riemannian mani- 
fold (M 2 (R),g). 

Note that such Riemannian structures defined by the 
metric tensor in a chosen representation are in fact 
closely related to classical thermodynamic fluctuation 
theory Q and critical phenomena. The probability dis- 
tribution of thermodynamic fluctuations in the equilib- 
rium intrinsic space naturally characterizes the invari- 
ant interval of the corresponding thermodynamic geom- 
etry, which in the Gaussian approximation reads Q(x) = 
A exp[—^gij(x)dx l (g> dxP], where the pre-factor A is the 
normalization constant of the Gaussian distribution. The 
associated inverse metric may easily be shown to be the 
second moment of the fluctuations or the pair correlation 
functions and thus may be given as g 1 ^ =< X l \X^ >, 
where {Xj}'s are the intensive thermodynamic variables 
conjugated to {x 1 }. Moreover, such Riemannian struc- 
tures may further be expressed in terms of any suit- 
able thermodynamic potential obtained by certain Leg- 
endre transform(s) which corresponds to general coordi- 
nate transformations on the equilibrium intrinsic mani- 
fold. Our geometric formulation thus tacitly involves a 
statistical basis in terms of a chosen ensemble, although 
the analysis has only been considered in the thermody- 
namic limit of QCD. 

We wish to mention that the thermodynamic curva- 
ture corresponds to the nature of the correlation present 
in the statistical system. This in particular will imply 
that the scalar curvature for two component systems can 
be thought of as the square of the correlation length at 
some given QCD transition temperature Tc(^i, ^2) to 
be R(fii,fi2) ~ S, 2 , where we may identify the £(Tc) to 
be the correlation length of the thermodynamic system. 
In turn, it is not surprising that the geometric analysis 
based on a Gaussian approximation to the classical fluc- 
tuation theory precisely yields the correlation volume of 
QCD. This strongly suggests that even in the context of 
chemical reactions in any closed systems, non-zero scalar 
curvature might provide useful information regarding in- 
teractions between various components of a QCD config- 
uration. 

Note further that the relation of a non-zero scalar cur- 
vature with an underlying interacting statistical system 
remains valid even for higher dimensional intrinsic man- 
ifolds and the connection of a divergent scalar curvature 
with phase transitions may accordingly be divulged from 



the Hessian matrix of the considered QCD free energy. 
It is worth to mention that our analysis takes account of 
the scales that are larger than the correlation length and 
considers that only few QCD microstates do not domi- 
nate the macroscopic quantities. In particular, we shall 
focus on the interpretation that the underlying QCD free 
energy includes contributions from a large number of mi- 
crostates and thus our description of the geometric QCD 
thermodynamics. With this general introduction to the 
thermodynamic geometry defined as the Hessian function 
of the free energy, let us now proceed to investigate the 
free energy of hot QCD and its thermo-geometric struc- 
tures. We now proceed to compute QCD free energies 
in the required form and thereafter study the thermody- 
namic geometry of the quark theories arising from the 
free energy of hot QCD. 



III. THE FREE ENERGY OF HOT QCD 

To obtain the free energy as a function temperature 
and quark chemical potentials, we consider the quasi- 
particle model of hot QCD. We consider two cases in the 
present paper viz. the 2-flavor QCD with two distinct 
chemical potentials (fi ~ (pi, {12)) an d the 3-flavor QCD 
with two distinct chemical potentials (fj,i corresponds 
to up and down quarks and [i 2 corresponds to strange 
quark). In order to obtain the expression for the free en- 
ergy we consider the following equilibrium distribution 
functions for gluons and quarks: 



f 



(exp(-/3e p )-l)' 

Jeq (exp(-/3e p - pp) + 1) (exp(-/3e p + pn) + 1) 

(5) 

For to obtain the free energy one utilizes the well 
known thermodynamic relation, 

F = -^EM^) + M^4 (6) 

p 

which is nothing but the negative of the pressure of hot 
QCD. Here Z g and Zm are the grandcanonical parti- 
tion functions corresponding to the pure gluonic part and 
matter sector respectively and have the following well 
known form in terms of the gluons and quark-antiquark 
distributions functions: 



HZ a ) 
ln(Z M ) 



- ln(l - exp(l - exp(-/3e p ))), 
ln(l + cxp(-/3e p — /i)) 
+ ln(l + exp(-/3e p +//)). 



(7) 



We now turn our attention to compute the free energy 
in the small chemical potential limit (j3fi << 1). We 
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shall determine the free energy only up to quartic terms 
in The expression thus obtained reads in compact 
notations 



F = T^ ao+Cl ^|+ C2 |i+ c l|i+^||)+T 4 0( M VT 6 ). 

We shall firstly consider the first three terms for our 
analysis which we denote as case (A) and the full expres- 
sion will be denoted as case (B). We shall study these two 
cases systematically and thereby show that the geometric 
quantities remain intact at zero chemical potentials. 

On comparison of the free energy displayed in Eq|8]for 

2- and 3-flavors with the temperature dependence of the 
corresponding pressure either in lattice QCD [lfl [ill [2(| 
or improved pQCD [H], HHU one may conclude 
that the coefficients of the and /i 2 should be negative 
in the corresponding pressure. We shall utilize this to fix 
the sign of the various terms in the Taylor expansion of 
the free energy as a function of the chemical potentials 
(/ii and /x 2 ). The various coefficients in Eq[5]for 2 and 

3- flavour QCD are as follows: 

for 2-flavor QCD, 



a 



87T 



1 



42 
32 



, ci = c 2 = 1/2, c\ 



An 2 



.(9) 



and for 3-flavour case the coefficient reads 
- 2 t 63 



a 



45 V 32 
1 

'An 2 ' 



,ci 



l,c 2 = 1/2,4 = - — 



1 

2n 2 ' 



(10) 



Note that, all the above coefficients are dimensionless. 
Let us now turn our attention to analyzing the free energy 
of hot QCD near Tc and trade off the Tc dependence in 
favour of quark chemical potentials following the relation 



(TP)2- 



(11) 



where &i = a 2 = —0.07 for the 2-flavor case and <ii = 
h 2 = -0.114 for 3-flavor QCD. We choose T c ° = 0.203 
Gev for N f = 2 QCD and T c ° = 0.197 Gev in 3-flavor 
QCD [26] for our analysis. After substituting for the 
transition temperature, Tc as a function of quark chem- 
ical potentials in the free energy, we thus obtain the fol- 
lowing polynomial form for the free energy: 

F(fj,x,fx 2 )- = a (b + ax^ 2 + a 2 ^l'f + (b + ai^l + a 2 /J,l) 2 

x(d^ + CBA»l)+ciV + 4 2) /4. (12) 

where b — T®, a\ — a%/b and a 2 = a 2 /6. Both a\ and a 2 
have the dimension of 6 _1 (GeV). The free energy will be 



measured in units of GeV 4 . Note that the components of 
the thermodynamic metric tensor are in units of GeV 2 , 
hence the thermodynamic curvature is in units of GeV 4 , 
and the chemical potentials \i\ and /i 2 will be in units of 
GeV . 

The above expression for the free energy is now only a 
function of two quark chemical potentials (^i,^ 2 ) since 
the temperature dependence has been traded off in favor 
of quark chemical potentials. The underlying covariant 
thermodynamic geometry arising from this free energy 
of hot QCD can be determined by employing the gen- 
eral formalism quoted in the previous section, which we 
shall do systematically in Section 4. Let us first discuss 
how one incorporates thermal fluctuations into the free 
energy and later on we shall study the contributions to 
the geometric structure of free energy. 

Let us now compute the contribution to the thermo- 
dynamic geometry from the thermal fluctuations. Let us 
consider case (A) for the 2-flavor QCD first. The free 
energy defined as the negative logarithm of the partition 
function takes the general polynomial form 



fW 2 ) :=«o + C i(f ) 2 +c 2 (f ) 2 - 



(13) 



Since the specific heat as the thermodynamic param- 
eters at constant volume may be defined as the second 
partial derivative of the free energy with respect to the 
temperature, this enforces that the associated leading or- 
der free energy with the inclusion of the thermal fluctu- 
ations after substituting for Tc as a function of chemical 
potentials reads 



F (^1,1^2) 



a 



cinl + C 2^2 



ainl - 



2\2 



2 K {b + ai ^\ 



■ C2M2 



a 2 ^)2 



)• (14) 



Let us turn our attention to explore the thermody- 
namic geometry of the equilibrium configuration space 
of the hot QCD, arising from the free energy expression 
at the various order perturbative contributions as well as 
that of the thermal fluctuations at the leading order free 
energy. In order to analyze the geometric nature of the 
hot QCD, we shall focus our attention on the realistic 
values characterizing certain specific phases of the QCD. 
In particular, the examples of interest are the case of two 
flavor QCD or that of 3-flavor QCD. As the exact expres- 
sion for the components of the metric tensors, associated 
determinants and that of the scalar curvatures are rather 
involved with general coefficients, thus we shall provide 
the geometric expressions only for the realistic parame- 
ters of the free energy which describes the present nature 
in hot QCD near the Tc- For simplicity, we have not 
explicitly considered the mass of the quarks but on gen- 
eral ground, our geometric conclusions may be expected 
to retain the same conclusions. Moreover, this is an in- 
teresting path to pursue and, for instance by arguments 
similar to the one given above we can conclude that the 
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general results with the inclusion of the thermal fluctua- 
tions agree with the results for arbitrary perturbatively 
corrected free energy in the framework of the theories of 
the quasi-particles. We shall now proceed to study the 
thermodynamic geometry in Case(A) and Case(B) one 
by one in the 2- and 3-flavor QCD. We shall henceforth 
suppress the units of the thermodynamic metric (GeV 2 ) 
and thermodynamic curvature (GeV 4 ). 

IV. 2-FLAVOR QCD 

In this section, we shall present the essential features of 
thermodynamic geometries and thereby put them into ef- 
fect for the 2-flavor hot QCD. In particular, we focus our 
attention on the geometric nature of QCD in the neigh- 
borhood of Tc with given quark chemical potentials in- 
troduced in the framework of quasi-particlc theories. As 
stated earlier the thermodynamic metric in the chemical 
potential space is given by the Hessian matrix of the free 
energy with respect to the intensive variables which in 
this case are the two distinct chemical potentials carried 
by the quarks. 

A. Case(A) 

Considering EqfTJlin the neighborhood of Tc and sub- 
stituting for Tc as function of [i\ and /12, one obtains the 
following expression: 

F(/ti,/i 2 ) : = a (b + aifil + a 2 fJ,l) 4 + {b + ami + a 2 nl) 2 
x(ci/4 +c 2 /4)- (15) 

To obtain the thermodynamic metric tensor in the 
chemical potential space, we employ the formula in Eq. 
[3J which leads to the following polynomial expression for 
the components of the metric tensor: 



Xij (Mi > M2) = Xij (^2, Mi); where i,j are the flavor indices 
associated with the 2-flavor QCD configuration, see for a 
phenomenological example [ij. Furthermore, the deter- 
minant of the metric tensor turns out to be polynomial 

ff(/*i.J*a) - -5.81(^° + M2°)-0.39(m?+m!)+0.02 
+37.85(MiM| + M? + i4) - 29 M(i4l4 
+/^ + M!)-58.07(/x?/4 + Mi + M2) 
+56.78/4/4 + 26.11/4/4 " 21.22(/4/4 
+/4 + m!) + 7.83(mi/4°+Mi°+m!) 

+i.3i(/4 2 + /4 2 )-7.07(/4 + /4) 

+9.460*5 +/4)+2.53(/4 + /4) 
+19.59(/4/4 + /4/4). (17) 



Thermodynamic metric in 2-flavor QCD 
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FIG. 1: The determinant of the thermodynamic metric in 2- 
flavor QCD in the chemical potentials surface. Note that we 
measure the chemical potential in GeV. 
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1.12/4 + 5.60/4 + 6-72/4/4 - 6.48/4/4 
-3.89/4/4 - 2.21/4 - 0.74/4 + 0.13 
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-1.48/ii/i 2 + 4.48/4/i 2 + 4.48/ii/4 
-2.59/4/i 2 - 2.59(iif4 ~ 5 - 18 /4/4 
1.12/4 + 5.60/4 + 6.72/4/4 - 6.48/4/4 



-3.89/4/4 



2.21/4 



0.74/4 



0.13 



-3.02/4 - 0.43/4- 



(16) 



We thus see that the components of metric functions 
indeed satisfy (i) (fii, Ma) = 9n a ,fn (M2, Mi) and (h) 

9 mi, M2 (Mi = M2) = 5^1,^2(^2,^1). Our thermodynamic ge- 
ometry predictions thus obtained indicate that the nature 
of two point correlation functions of such 2-flavor QCD 
must satisfy (i) and (ii). In particular, this determines 
the nature of associated quark number that the suscepti- 
bility tensor satisfies: Xijif^i, M2) = Xji(^2, Mi) = X an d 



Employing the formula displayed in Eqf4l we obtain 
the thermodynamic curvature as 



fl(Mi,M2) = -^{0.31(m1°+M2°) + 1-37(Mim1 + Mi+M2) 

g 2 

+0.82(m1m2° + + /4) + 0.27(m1 2 mI + Mi 

+/4 2 ) + o.oo4(/4 + (4) - 1.18(^(4 + /4 

+M2) + 1-55(m 2 mI + Mi + l4) 

+z.io{i4(4 + (4 + &) -om^l 

-1.77(4(4 + 0.46(/4/4 + /it + M2) 

-i.02(/4/4° + /4° + /* 2 ) - o.i7(/4 2 + /4 2 ) 

-3.42m?m| + 0.15(/4 + M2) - 0.29(/4 + Ma) 
-0.04(/4 + /4) - 2.56(m^ + M? + M2) 
+0.04(/4 4 + M 2 4 )}- 

(18) 
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FIG. 2: Thermodynamic curvature in 2-flavor QCD in the 
chemical potentials surface. Note that we measure the chem- 
ical potential in GeV. 



The determinants of the thermodynamic metric and 
curvature in the chemical potential space are shown in 
Figs 1 and 2 respectively. We denote (/xi,/^) = (u, v), 
and henceforth focus our attention in the physically in- 
teresting domain near Tc, where both u and v lie be- 
tween 0.0 to 0.2. In this case, we see that there are 
two distinct curves of certain bumps of varying lengths 
whose size depends on the domain chosen in the (u, v) 
space. The significance of these bumps is that they show 
thermodynamical interactions in the underlying system. 
For the determinant of the metric tensor, when plotted 
against the chemical potentials, we find that it is a reg- 
ular increasing function of the chemical potentials in the 
considered physical domains of hot QCD. 



the thermodynamic geometry remains conserved as that 
in the Case(A). This in fact remains true in general for 
any symmetric function F ((11,(12) whose Hessian matrix 
defines the metric tensor of the underlying intrinsic 
Riemannian manifold. The determinant of the metric 
tensor turns out to be a polynomial function in the 
chemical potentials {/!i,/!2}, which in turn is given by 

g(M,H2) = 19.59(^ + ^)-5.81(^ n + M2°) 
-7.42 Gu? + M f)-0.43( M f + /4) 

-29.03(/i?/xl + M1M2) + 2.75(/4 + /4) 
+0.02 - 58.07(/4mI + Mi/4) 
+7.83(/x^ + /i 1 ^) + 26.12/i^ 

+6.50(4(4 + 9.60(^1 + mI) 
-24.97(^2 + Mi/4) + l-31(Mi 2 + Ma 2 ) 
+39.96(/z 2 ^ + M1M2) + 60.73m?M2- ( 20 ) 
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B. Case(B) 

Considering the expansion of the free energy given in 
Eq. 1121 and employing Eq. [3 for non vanishing c\,cl 
defined in Eq. [51 the following expression for the compo- 
nents of the metric tensor may be obtained: 

9^,^ = 6. 72(4(4 - 2. 52(4 - 0.74/i 2 + 0.13 + 5.60/4 
+1.12(4 - 3.02/x? - 0A3f4 - 6.48m1mI 
-3.89 M 2 M2, 

9 nuv* = -1.48/xxM2 + 4.48/x?M2 + 4.48/ximI - 2.59MiM2 
-2.59/xi/i^ - 5.18miM2) 
= 6.72f4f4- 2.52(4 -0.74//? + 0.13 + 5.60/4 
+1.12/4 - 3.02/^ - 0.43juf - 6.48/4/4 
-3.89/4mI- (!9) 

We see further that in this Case(B), the 
symmetry(/ii — > — Mi) an< i exchange symmetry of 



FIG. 3: The determinant of the thermodynamic metric in 2- 
flavor QCD as a surface plotted against the chemical poten- 
tials. Note that we measure the chemical potential in GeV. 



We see that the determinant of the metric tensor re- 
mains non-zero in the small chemical potential limit and 
thus defines a non-degenerate thermodynamic geometry 
near Tc- Finally, we may easily obtain the underlying 
thermodynamic curvature, and it reads as follows: 

R((H,(i 2 ) = -^(Vo4(/4 4 + /4 4 )-3.09(/4 M 4 + / l2 V4) 

+0.28(^ + M2°)+0.15( Ml + M2)) 
+0.01(/i 2 + (4)) + 1.87(/x?Mi + M2M1) 
-0.04(/4 + (if) + 4.21(m?M2 + M2M1) 
-1.13(M?Mi° + M2M1 ) + 0.27(m 2 m 2 2 
+/4V 2 ) + 0.82(/4/4° + /4Vi) 
-4.23/i?M2 + l-37(MiM2 + mImi) 
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-1.0/4/4 - 0.27(/4 



+/4/4) - 0.16(/4 2 + /4 2 



+l4t4) - 2.49/4/4 



l4) + 0.59(/4/4 
l-51(/i?M2 



(21) 
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FIG. 4: Thermodynamic curvature in 2-flavor QCD in the 
chemical potentials surface. Note that we measure the chem- 
ical potentials in GeV. 

The behavior of g and R as a function quark chem- 
ical potential potentials in Case (B) for 2-flavor QCD 
are shown in Figs. 3 and 4. Here the conclusions to be 
drawn in the range of the chemical potentials considered 
in the previous case remain the same, except for the fact 
that the curves of the interactions are aligned with an en- 
hanced magnitude, while the shape of the determinant of 
the metric tensor plotted against the chemical potentials 
remains (exactly) the same as in the previous case. 

We see however that the components of the metric ten- 
sor at zero chemical potentials take the values = 
0.13 = 5p 2 .p 2 and g Ml , M2 = 0, which effectively describe 
nothing more than the diagonal and, respectively, off di- 
agonal quark susceptibility tensors. We find further that 
the determinant of the metric tensor reduces to the value 
of g = 0.02 at zero chemical potentials, while the scalar 
curvature becomes zero at this point and thus the un- 
derlying statistical system at this point becomes non- 
interacting. We further see that these diagonal and off 
diagonal quark susceptibility tensors, as well as the de- 
terminant of the metric tensor and the scalar curvature 
at zero chemical potentials, remain the same under the 
(^^-contributions. 

In the case of the 2-flavor QCD free energy, we note 
that the degree of the determinant of the thermodynamic 
geometry in the chemical potential space remains intact 
for the Case(B) as that of the Case(A). We further see 
that the associated thermodynamic curvature never di- 
verges for any values of the chemical potentials, and tends 
to a vanishingly small quantity when either one of the 
chemical potential corresponds to a large value. Our pre- 
dictions can thus easily be generalized to the case when 



the free energy may be expressed as a finite polynomial 
expression of two chemical potentials. We have further 
noticed that the thermodynamic geometry of the free en- 
ergy preserves all the symmetries of the free energy. Let 
us now proceed to study the effects of thermal fluctu- 
ations to the thermodynamic geometry of case (A) in 
2-flavor QCD. 



C. Thermal fluctuations 

We will now discuss the thermodynamic geometry of 
hot QCD with the inclusion of thermal fluctuations about 
an equilibrium. As it is well known, any thermodynam- 
ical system considered as a canonical ensemble has loga- 
rithmic and polynomial contributions to the free energy 
These considerations apply as well to quasi-particles 
models of QCD, and the specific forms of the logarith- 
mic and polynomial contributions may easily be calcu- 
lated for a wide class of particle theories. In order to 
the analyze essential geometric nature of hot QCD under 
thermal contributions, we shall consider the case of the 
quadratic corrected free energy. Further it is straight- 
forward to see that the case of quartic corrected free en- 
ergy in this framework leaves the qualitative nature of 
results untouched. This approximation is valid only in 
the regime where thermal fluctuations are much larger 
than quantum fluctuations. 

With the above consideration in the leading order 
quadratic approximation, the corresponding free energy 
of two flavor QCD reads 



do 



Ci/4 



C2M2 



(6 + 

-~ln(— 

2 \b 



ci/4 



a 2 /4) 2 
-c 2 i4 



(22) 



Employing the formula to determine the metric ten- 
sor and thermodynamic curvature displayed in Eqs[H the 
components are the metric tensor in the scale of 10 10 
come out to be 



9n 



{-0.17(/4 - /4) - 1.33(/4° 



/4°) 

+9.96/4 - 36.26/4 - 3.55/4 - 5.82/4 
+2.67(^/4 - /4/4) - 102.80(/4/4 - 0.5/4/4) 
-78.93/4/4 - 36.27/t? - 51.40/4 - 5.82/4 
-1.33/4° - 6.38/4/4 ~ 98.86/4/4 

-4.oo(/4/4 - /4/4) - 9.37/4/4} 

x (-203 + 340/4 + 340/4r 4 (/4 + /4)~ 2 \ 
Wa = 2/ 11 // 2 {1.13(/4 + /4) + 23.12(/4 + /4) 
+1.33(/4 + /4) + 51.40(0.50/4 + /4) 
+69.36(/4/4 + /4/4)8.01(/4 + /4) 

+5.34(/4/4 + /4/4) + 51.40/4/4 - 0.17} 
x (-203 + 340/4 + 340/4r 4 (/4 + vl)^ \ 
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= {-0.17(^ - (i\) - 1.33(A»i° - a4°) 
+9.96/4 " 36.26/4 - 3.55yu| - 5.82/4 
+2.67(/4/4 " a4/4) - 102.80(m^ 
-0.5/4(4) - 78.93/4/4 ~ 36.27/4 
-51.40(4 - 5.82(4 - 1-33(4° - 6.38/4/4 
-98.86/4/4 - 4.00(/x?m1 - /4m?) 
-9.37/4/4} x (/4 + M2)" 2 
(-203 + 340/i? + 340/4)~ 4 - (23) 

It is thus evident that the symmetry of the metric ten- 
sor remains preserved under thermal fluctuations. In this 
case, after incorporating the Log contributions to the free 
energy in the case of 2-fiavor QCD, we obtain the follow- 
ing positive definite expressions for the determinant of 
the thermodynamic metric: 

gipuin) = lO- 3 (Vo4(/i 2 -Mi) + 34.38(/4 + /4 ) 

+29.64(/4 + (4) + 58.41(/i? + (4) 
+3.69(/4 + (4) + 27.7304V? + MiVI) 
+343.88 (/zVi + /4/4) + 88.92(/4/4 
+(i\(4) + 350.50(4(4 + 4.62(//J 2 + /4 2 ) 
+0.13(/i! 4 + m") + 69.32(/4m! + M1M2) 
+2.75(/4Vi + a4/4°) + 171.94(^/4 
+/4/4) + 92.43/^/4 + 233.67(/4/4 
+/4/4) + 4.59( M Vf + /4/4) ~ 0.035 

+om{(4(xf + /4V?) + 7.38MV2) 

x K0.5(/4 + /4))~ 2 (0.20 - 0.34/4 
-0.34/4)- 7 ). (24) 

We see in this case again that the thermodynamic cur- 
vature turns out in the form of the ratio of two polyno- 
mial expressions, which may be written as 

R((ii,(i 2 ) = -4(0-203 - 0.34(/4 + /i 2 ) 
9 

24 

a i,jf4^2> ( 25 ) 

{i,j=0|i+j<24} 

where the coefficients {a^} of the polynomial expres- 
sion appearing in the numerator may easily be tracked 
from the corresponding points on the figure 6. In this 
case, we should note that the determinant of the metric 
tensor plotted against the chemical potentials becomes a 
smoother function in comparison to that of the case with- 
out logarithmic correction. However, we observe that it 
reaches a magnitude of 10 9 near the (u, v) = (0.2,0.2). 
Further, we find that the domain of the interaction has 
shifted towards the origin. In particular, we see that such 
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FIG. 5: The determinant of the metric tensor in 2-flavor QCD 
under thermal corrections in the chemical potentials surface. 
Note that we measure the chemical potential in GeV. 



interactions are only present in the chemical potentials 
range 0.1 < (ii < 0.2. 
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FIG. 6: Thermodynamic curvature in 2-flavor QCD under 
thermal corrections in the chemical potentials surface. Note 
that we measure the chemical potential in GeV. 

Finally, we wish to point out a general feature of our 
geometric considerations, i.e. that any polynomial free 
energy after including the logarithmic contributions is 
going to generate certain bump(s) in the intrinsic Rie- 
mannian space of chemical potentials for a given non-zero 
Tq around the equilibrium QCD thermodynamic config- 
uration. Furthermore, following our general procedure 
it may explicitly be verified, in the same way as in the 
cases of quadratic or quartic contributions to the free en- 
ergy, that the logarithmic contributions to the free energy 
do not spoil the symmetry of the determinant and that 
of the associated scalar curvature of the intrinsic space 
spanned by the quasi-particle chemical potentials. Such 
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symmetries are indeed the artifact of the choice of the pa- 
rameters appearing in the corresponding free energy such 
that they correspond to the symmetric free energy in the 
chemical potentials. Due to the inherent symmetries in 
the definition of the quasi-particle free energies, it is pos- 
sible in principle that we can give up manifestly invariant 
expressions for the determinants and the scalar curva- 
tures. In particular, we find that for all quasi-particle 

free energy F(/ti,/t 2 ) := <H,jlAl4 witn a i,i = a j,^ ^ 
turns out that the determinant and associated scalar cur- 
vature of the intrinsic Riemannian manifold as a function 
of the quark chemical potentials are symmetric functions 

0(Mi>M2) = 9(1*2, Mi) and i?(/ii,M2) = #0"2,Mi) 01 tnc 
chemical potentials. 

V. 3-FLAVOR QCD 

As in the previous section, we shall analyze the free 
energy as a function of quark chemical potentials with 
or without thermal corrections, and thereby systemat- 
ically explore the covariant metric tensor, determinant 
and scalar curvature of the thermodynamic geometry of 
3-flavor QCD. 

A. Case(A) 

In what follows next in this subsection, we present the 
geometric analysis for the realistic quasi-particle model 
of QCD near the Tc- The thermodynamic geometry of 
interest may thus easily be obtained, as before, from the 
Hessian matrix of the free energy with respect to the 
quark chemical potentials. In this case the components 
of the metric tensor reduce to 

gp ulil = -7.59//? - 2.30(4 + 0-26 + 33.76/4 
+6.08/4 " 32.50/4 + 38.51/4/4 
-69.65/4/4 - 41.79/4/4 ~ 4.64/4, 

fa,*, = -4.60/xi/* 2 + 25.67/4^2 + 24.33/i X /4 
-27.86/4^2 - 55.72/4/4 ~ 27.86/ti/t^ 
= -6.22/4 - 2.30/4 + 0.22 + 6.41/4 

+28.75/4 - 4-64/4 + 36.50/4/4 
-41.79(4(4 - 69.65/4/4 - 32.50/4 (26) 

We see further that the components of the metric ten- 
sor at zero chemical potentials take the values g^,^ = 
0.26, <7 M2j M2 = 0.22 and g^.^ = 0, which in effect re- 
spectively describe the diagonal and off diagonal quark 
susceptibility tensors. We further find that the determi- 
nant of the metric tensor reduces to the value of g = 0.06 
at zero chemical potentials, while the scalar curvature 
vanishes at this point and thus the underlying statisti- 
cal system approaches a non-interacting system. On the 
other hand, the determinant of the metric for this 3 flavor 
case is as follows: 



</(/ii,M2) = 2263.60(/4/4 + /4/4) + 905.44(/4/4° 

+/4/4°) + 3018.13/4/4 - 2.30(/4 + (4) 
+0.06 + 49.23/4/4 ~ 380.47(/4/4 + /4/4) 
+1157.67/4/4 + 1254.19/4/4 + 150.91(/tJ 2 
+/4 2 ) - 1792.58/4/4 - 1690.20/4/4 
-3516.91/4/4 ~ 3448.66/4/4 + 26.65/4 
+23.19/4 - 113.54/4 - 134.84m? 
+278.58/4 + 326.69/4 + 1806.89/4/4 
-331.22/4° - 365.34/4°. (27) 
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FIG. 7: The determinant of the thermodynamic metric in 3- 
flavor QCD in the chemical potentials surface. Note that we 
measure the chemical potential in GeV. 



In this case we see that g((i\,(i2) 7^ <?(/i2, Mi), which 
follows from the fact that the 3-flavor QCD free energy 
(displayed in Eq|12p is not symmetric under the exchange 
in the chemical potentials. The same remains true for the 
associated curvature, which may be seen from 

#(//!, /t 2 ) = -^-910.89/4/4-374.91/4/4° 

-1182.43/4/4 + 190.37/4/4 2 
+554.91/4/4° + 0.078/4 + 0-08/4 
+897.84/4/4 + 163.36/4 2 /4 - 2.58/4/4 
+25.72/4/4 - 338.28/4°/4 + 870.84/4/4 
-122.29/4/4 + 506.30/4°/4 + 26.01/4/4 
+282.79/4/4 - 865.09/4/4 ~ 55. 24//} 2 
+574.23/4/4 - 64.40/4 2 - 1.28/4 
+585.53/4/4 + 8-82/4 + 22.57/4 4 
-1.31/4 + 27.97/4 4 + 8.54/4 
-119.06/4/4 + 299.74/4/4 
-31.25/4 - 29.65/4 - 180.44/4/4 
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+61.61/4° + 55.96/4° . (28) 
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FIG. 8: Thermodynamic curvature in 3-flavor QCD in the 
chemical potentials surface. Note that we measure the chem- 
ical potential in GeV. 

From Figs. 7 and 8, we see that the determinant of 
the metric tensor plotted against the chemical potentials 
shows two lines of minima and the maximum lying in 
between the minima occurs when either one of the chem- 
ical potentials reaches a non-zero constant value, while 
the other is fixed at the origin. In this case we further 
see that there are two distinct bumps of varying heights 
instead of the curves of interactions. The height of the 
bumps depends on the domain chosen in the (u, v) space. 
The significance of these bumps is that they show ther- 
modynamical interactions in the system of the 3-flavor 
QCD. 

B. Case(B) 

In this subsection, we shall investigate the thermody- 
namic geometry corresponding to the free energy given 
by EqEl It turns out here that the components of the 
metric tensor thus obtained are as follows: 



remains the same, except the fact that the number of 
bumps is increased, which thus indicates stronger inter- 
actions in the internal space of the quark chemical po- 
tentials, while the shape of the determinant of the met- 
ric tensor plotted against the chemical potentials is a bit 
flatter near the (u, v) — (0.2, 0.2) than the previous case. 
The determinant of the metric tensor turns out to be 
polynomial 

5 (/ii,/i 2 ) = -331.21/4° - 365. 34/4°- 138.74/4 
-115.38/4 - 2.22/4 - 2.42/4 + 0.06 
-1690.20/4/4 - 1792.58/4/4 + 23.89/4 
-3448.66/4/4 + 28.06/4 - 3516.91/4/4 
+905.44/4/4° + 2263.60/4/4 + 150.91/t} 2 
+2263.60/4/4 + 279.99/4 + 905.44/4°/4 
-412.95/4/4 + 329.52/4 + 1289.51/4/4 
+1870.46/4/4 + H90.46/4/4 + 55.71/4/4 
-388.38/4/4 + 3018. 13/4/4 
+150.91/4 2 . (30) 
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FIG. 9: The determinant of the metric tensor in 3-flavor QCD 
as a function of the chemical potentials. Note that we measure 
the chemical potential in GeV. 



gn ulll = -38.51/4/4 + 6.08/4 ~ 69.65/4/4 Note f urtner in tms case > we a S ain see tnat 9(^1,^2) 

—8 20 2 — 2 30 2 + 61 + 33 76 4 ff(M2>Mi)) which follows directly from the fact that 

1 2 4 ' 2 g ' 1 the underlying free energy with (^^-contributions (dis- 

— 41.79/i x /t 2 — 32.50/ij — 4.64/x 2 , played in Eq[8]) remains non-symmetric under the ex- 

5mi,M2 = — 4.60/ii/i2 + 25.67/tJ/ii + 24.33/ii/t2 change of the quarks chemical potentials. Thus the same 

—27 86/4/t 2 — 55 72/4/4 — 27 86/j,i/4 remains true for the associated curvature as the polyno- 

22 4 4 2 mial ratio of the Riemann tensor and the determinant of 

= -36.51/i x /x 2 + 28.75/i 2 - 41.79/z 1 /t 2 thc mctric tensor, which turns out to be given by 



-2.30/4 - 6.53/4 + 0.22 + 3.42/4 

-69.65/4/4 - 4 -64/4 - 32.50/4- (29) 4 

^2 



R(jii,[Mi) = — ji 59 - 73 ^2° + 8.32/4 + 8.71/4 + 0.09/4 
In this case the conclusion to be drawn in the range of 9 
the chemical potentials considered in the previous case +0.09/4329.65/4/4 684.70/4/4 — 1-34/4 
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-1.33/4 + 671.35/^/4 + 308.57/zfyu| 
-395.82/4/4° + 190 .37 M 2 M i2 _ i 04.36/4/4 

-1322.64/4/4 + 554.91/4/4° + 506.30(4° (4 
+897.84/4/4 + 870.83/4/4 - 952.42(4(4 
-30.34/4 - 27.68/4 ~ 354.26/4VI 
+28.95/4/4 - 133.49/4/4 - 212.79/4/4 
-137.42/4/4 + 29.18/4/4 - 2.82(4(4 
+22.S6/4 4 + 27.97/4* - 63.17/4 2 
-52.78/i} 2 + 163.36/4 2 /4 + 51.99/4°}. (31) 
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FIG. 10: Thermodynamic curvature in 2-flavor QCD in the 
chemical potentials surface in case B. Note that we measure 
the chemical potential in GeV. 

In these cases when the free energy may be treated 
as a function of the chemical potentials, we see that the 
thermodynamic space spanned by the chemical potentials 
remains completely regular except at those points for 
which the determinant of the metric tensor vanishes on 
the intrinsic manifold. This intrinsic space of the chemi- 
cal potentials thus turns out to be a well-defined, stable 
configuration for all {(11,(12} such that the the determi- 
nant of the associated metric tensor remains positive def- 
inite in the domain of interest. We see in this case that 
the underlying diagonal quark susceptibility tensor does 
not remain the same at zero chemical potentials, and 
in particular, we find that = 0.61, = 0.22 

and g^,^ = while the determinants of the metric ten- 
sor and the scalar curvature remain the same under the 
(^^-contributions at zero chemical potentials. 

C. Thermal Fluctuations 

We will now discuss the thermodynamic properties of 
3-flavor QCD under thermal fluctuations, treating the 
system as a well-defined ensemble. We allow for small 
thermal fluctuations in the system considered, and study 
the thermodynamic geometry of 3-flavor QCD in the lines 



with our treatment of the usual two flavor QCD described 
earlier in the previous section. For simplicity we shall 
work in the leading order approximation with the cor- 
responding free energy as that of two flavor QCD, with 
the understanding that the parameters appearing in the 
free energy are given as that of the quadratic case of 3- 
flavor QCD. A straightforward computation yields that 
the components of the metric tensor in the scale of 10 10 
read 

Wi = -2(-197 + 578/4 + 578/ l2 T 4 (2/4 + /4)- 2 
{7.50/4 + 3.80/4 + 0.30/4 ~ 0.15/4 
-2.56/4 + 182.75/4/4 - 28.69/4/4 
+78.13/4/4 - 22.32/4° - 416.12/4 
-343.63/4 - 11.16/4/4 ~ 91.30/4/4 
-19.06/4/4 - 274.74/4/4 ™ 20.83/4 
+62.18/4/4 + 100.45/4/4 + 33.48/4/4}, 

Wl = 4/ ilM 2(-197 + 578 Ml + 578/4)- 4 (2/ t 2 + / t 2)- 2 
{7.61/4 + 1-77/4 + L77 /4 + 27.68/4 
+78.13/4/4 + 107.83/4/4 + 33.48/4 
+385.68/4 + 89.29/4/4 + 126.27/4/4 
+385.68/4/4 + 143.05/4 + 22.32/4/4 
-0.15}, 

W2 = (-197 + 578/4 + 578/4)- 4 (2/4 + /4)- 2 
{80.72/4 + 107.83/4 + 0.30/4 ~ 0.15/4 
-11.16/4° + 6.53/4 + 1911-89/4/4 
+368.99/4/4 + 66.97/4/4 ~ 69.97/4° 
-340.04/4 + 75.53/4 - 100.54/4/4 
+1278.78/4/4 + 17.29/4/4 + 437.31/4/4 
+ 15.52/4 + 801.80/4/4 + 178.58/4/4 
+89.29/4/4}- (32) 

The determinant of the metric tensor turns out to be 
given by the following polynomial expression: 

g((i 1: (i 2 ) = {0.067/4 + 0.109/4 + 0.00008/4 

+0.0001/4 + 0.0037/4 2 + 0.044/4° 
+0.008/4 + 0.032/4/4 2 + 0.194/4/4° 
+0.485/4°/4 + 3.292/4/4 + 0.409/4/4 
+12.120(4(4 + 3.786/4° + 3.160/4 
+6.190(4(4 + 0.529(4 + 0.032/4 4 
+0.778/4 2 + 11.250/4/4 + 0.394/4/4° 
+3.490/4°/4 + 5.207/4/4 + 0.194/4 2 /4 
+0.030/4/4 + 0.816/4/4 + 0.026/4 
+0.485/4/4 + 5.404/4/4 + 0.647/4/4 
+0.959/4/4 + 5.680/4/4 + 2.319/4/4 
+1.068/4/4 - 5.76 x 10~ 6 } x 

^(0.197 - 0.578(/4 + /4)r 7 
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chemical potentials surface when the thermal contributions 



are included. Note that we measure the chemical potential in 

FIG. 11: The determinant of the metric tensor in 3-flavor GeV. 
QCD in the chemical potentials surface when the thermal con- 
tributions are included. Note that we measure the chemical 

potential in GeV. VI. DISCUSSION AND CONCLUSION 



In this case, we see that the determinant of the met- 
ric tensor plotted against the chemical potentials ac- 
quires two finite size bumps which were absent in the 
case without thermal corrections. Again the height of 
the bumps depends on the location being chosen. Fur- 
thermore, we find that the domain of the interaction has 
shifted towards the origin, as in the previously treated 
case of 2-flavor QCD. In particular, we see that there are 
two bumps of interactions present near the value of ei- 
ther chemical potentials [ii = 0.1. However, we observe 
that the strength of the interaction depends on the ex- 
act location chosen in the internal space, for example it 
may become as high as 10 4 or as small as 10 2 near the 

(U,V)= (0.1,0.1). 

In this case again it turns out that the thermodynamic 
curvature may be written in the same form as the ratio 
of two polynomial expressions, and in particular we find 
that this curvature is given by 

fl(Mi,M2) = -4(0-197 - 0.578(^ 2 + /i2)) 
9 

24 

&i,i/4/4 (34) 

{i,J=0|t+i<24} 

where the associated coefficients {bij} of the Riemann 
tensor appearing in the numerator of the scalar curva- 
ture may easily be read out from the corresponding Fig. 
12. Thus, in this section once again one can always find 
from the general considerations of our intrinsic Rieman- 
nian geometry as in the case of 2-flavor QCD that for 
any given polynomial free energy after including the log- 
arithmic contributions there exist certain bump(s) in the 
associated intrinsic Reimannian space. 



In this paper we have introduced an intrinsic geomet- 
ric notion to QCD thermodynamics, and thereby it has 
been applied to study the behavior of 2-flavor and 3- 
flavor QCD in several simple thermodynamic considera- 
tions. In particular, the free energy of the 2-flavor, as well 
as that of the 3-flavor QCD, with or without logarithmic 
contributions under thermal fluctuations, indicate an in- 
triguing relationships between the geometrical concepts 
of QCD thermodynamics, and in particular, the scalar 
curvature of the underlying metric tensor to the correla- 
tion volume. Thus the physical concepts of a phase tran- 
sition^), if any, as well as the geometric apprehension of 
the quasi-particles, non-ideal interactions are divulged. 
For two-component QCD as a thermodynamic system 
with two distinct quark chemical potentials, the associ- 
ated quarks susceptibility tensors may consequently be 
defined as the components of the thermodynamic met- 
ric tensor arising from the Hessian of the hot QCD free 
energy and thus represent the scalar curvature of the un- 
derlying geometry from the viewpoint of the statistical 
theory of quark states. We have explicitly determined 
such notions for a few simple cases, in particular for the 2- 
flavor and the 3-flavor theories, and in turn find that they 
indicate certain physical behaviors in the lieu of "macro- 
scopic/microscopic duality" relation(s), which are indeed 
in close connection with the "AdS/CFT correspondence" 
via an associated fcrmion number density of the under- 
lying effective field theory. In fact, this study provides 
convincing examples that our geometrical approach ana- 
lyzes the thermodynamics of hot QCD systems from the 
perspective of microscopic theory and thus can be applied 
to actually divulge important physical and chemical be- 
haviors of the quasi-particles. 

We have shown that the covariant intrinsic geomet- 
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ric notion about an equilibrium thermodynamics config- 
uration, when applied to the 2-flavor and 3-flavor quasi- 
particle QCD models considered as thermodynamic sys- 
tems, provided that they are an integral part of the en- 
vironment with which they are treated to be in equilib- 
rium, in the framework of an underlying ensemble the- 
ory, indicates in both cases a well-defined, curved, regu- 
lar intrinsic Riemannian manifold. Hence, the study of 
the equilibrium thermodynamic geometry as an intrinsic 
manifold arising from quasi-particles and the associated 
scalar curvature have led us to interesting insights into 
phase transitions and critical phenomena, if any, for such 
quasi-particles. However, it may be mentioned that a 
thermodynamic fluctuation approach is not viable for a 
most general non-perturbatively completed QCD in the 
framework of non-equilibrium statistical mechanics with 
an infinite extensive environment, which we would like 
to investigate in the near future. The basic idea involved 
here is that the positivity of the thermodynamic metric 
tensor cannot be ensured for the full range of all param- 
eters {/Ltj, T, V, S}, which is necessary for a local mini- 
mum of the energy for any given equilibrium configura- 
tion. Thus, it is plausible to investigate such notions in a 
restricted fashion where fluctuations of one or more ther- 
modynamic parameter is either frozen out or assumed to 
be negligible, see for example (2?| in the context of QCD. 

In this paper, we have introduced such geometries to 
the thermodynamics of hot QCD, and focused our atten- 
tion on the intrinsic Riemannian geometry in the neigh- 
borhood of Tc(/ii, 1x2) of the 2-flavor, as well as of the 
3-flavor hot QCD. The geometry thus defined on a two di- 
mensional manifold is solely spanned by the chemical po- 
tentials of quasi-particles and in turn determines the two 
point correlation functions and an associated correlation 
volume, near the considered transition temperature of 
the underlying microscopic theory. Note further that all 
those QCD-diagrams which just have one gluon exchange 
do not appear in our thermodynamic geometry, which is 
due to the fact that such diagrams vanish under color 
neutrality transformations. Whereas the corresponding 
diagrams with two gluon exchange remain non-zero, be- 
cause the fermion loops are even functions of the chemical 
potentials and thus contribute to the off-diagonal suscep- 
tibility tensors only at zero chemical potentials. We thus 
find that components of the off-diagonal susceptibility 
tensor macroscopically correspond to those of the ther- 
modynamic metric tensor {<?ij |i 7^ j}, which may further 
be envisaged in terms of the fermion number density as 
gij = Xij — §t^S O ur geometric study may thus de- 
scribes the mixing between different flavors induced by 
the resummation of quark loops along the soft, internal, 
gluon lines to certain orders of given QCD diagrams. 

It is worth to point out in the case of 2-flavor QCD 
that the components of the diagonal quark susceptibility 
tensor and those of the off-diagonal quark susceptibility 
tensor may easily be read off just from the components 
of the associated thermodynamic metric tensors at zero 
quark chemical potentials. In particular, we see that the 



diagonal quark susceptibility tensor components <7^ <)([tj 
take identical values, while those of the off-diagonal quark 
susceptibility tensor identically vanish. It is straightfor- 
ward to see that the determinant of the metric tensor 
at zero chemical potentials takes a positive definite value 
and thus defines a well-defined thermodynamic geome- 
try, while the associated scalar curvature at this value of 
the chemical potentials turns out to be zero, and thus 
the underlying statistical system at this point becomes 
a non-interacting system. It turns out that both the di- 
agonal and off-diagonal quark susceptibility tensors, the 
determinant and the scalar curvature of such associated 
metric tensor at zero chemical potentials remain intact 
under the (^^-contributions. 

Moreover, we observe the same for the case of 3-flavor 
QCD, i.e. that at zero chemical potentials the compo- 
nents of the metric tensor i(tl< reduce to identical val- 
ues and in fact vanish corresponding to the case of g^^ , 
which thus respectively describe the diagonal and the off- 
diagonal quark susceptibility tensors. In this case we find 
that the determinant of the metric tensor acquires a non- 
zero positive definite value, while the associate scalar cur- 
vature turns out again to be zero, and thus the underlying 
QCD system turns out to be a non-interacting statistical 
system at zero quarks chemical potentials. It is impor- 
tant to note in this case that the diagonal quark suscep- 
tibility tensor components at zero chemical potentials do 
not find the same values under the (^) 4 -contributions 
but rather we observe g fJ , 1 ,p 11 — 0.61, <7^ 2i ^ 2 — 0.22 and 
9m 412 = 0- This further implies that the associated 
quark susceptibility matrix is non-trivial, and thus the 
different quarks acquire non-identical two point correla- 
tion functions near Tq. However, the determinant of the 
metric tensor or the geometric invariant objects such as 
the scalar curvature at zero quarks chemical potentials 
indeed remain the same under (^) -contributions. 

It is worth to note that our geometrical results of QCD, 
when plotted in a physically interesting, near TjS, range 
against the chemical potentials are in turn very illumi- 
nating. In particular, for the case of 2-flavor QCD we 
find that there are two distinct curves of bumps of vary- 
ing strength of thermodynamical interactions, whose size 
thus depends on the chosen domain in the chemical po- 
tential space of the system. Furthermore, the determi- 
nant of the metric tensor, when plotted against the chem- 
ical potentials in the considered domains of hot QCD, has 
been observed to be a regular increasing function of the 
chemical potentials. Interestingly, the conclusions to be 
drawn in this range of the chemical potentials remain 
the same under the (^) -contributions and it turns out 
that the shape of the determinant of the metric tensor 
against the chemical potentials remains exactly the same 
as that without such contributions. However, the curves 
of the interactions present get aligned with an enhanced 
magnitude compared to that without (y ) 4 -contributions 
for 2-flavor hot QCD. It is further important to note 
that the determinant of the metric tensor when plotted 
against the chemical potentials becomes a smoother func- 
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tion under the thermal contributions. Under the latter 
the domain of the interactions is observed to be shifted 
towards the origin, and in particular they are dominant 
in the range when both chemical potentials lie within 
0.1 < m < 0.2. 

On other hand, the case of 3-flavor QCD in turn shows 
that the determinant of the metric tensor, when plot- 
ted against the chemical potentials, has two lines of the 
minima, and the corresponding maximum exactly occurs 
when either one of the chemical potential reaches the 
maximum, while the other is held fixed at the origin. In 
this case, it is observed that there are just two distinct 
bumps of the thcrmodynamical interactions of varying 
heights, instead of the two curves of 2-flavor QCD. Since 
such thermodynamical interactions are local in nature, 
thus the height of the observed bumps depend on the do- 
main chosen in the underlying space of chemical potential 
which characterizes our QCD system. We however find in 
the 3-flavor hot QCD that the shape of the determinant 
of the metric tensor appears flatter, and that the number 
of bumps is increased under (^) 4 -contributions, which in 
fact indicates that the underlying configuration charac- 
terized as an internal manifold of the chemical potentials 
is a strongly interacting statistical system. It turns out 
that the determinant of the metric tensor, with the in- 
clusion of thermal fluctuations, acquires two finite size 
bumps when plotted against the chemical potentials. In 
particular, we observe that there are two bumps of inter- 
actions present when either one of the chemical poten- 
tials reaches /x, = 0.1, whereas such bumps do not ap- 
pear in the determinant of the metric tensor in this range 
of chemical potentials. It may further be observed that 
the strength of the local thermal interactions becomes di- 
verse with the amplitude of 10 2 GeV 4 to l0 4 GeV 4 of the 
interactions near the (^1,^2) = (0.1, 0.1)GeV. 

We have thus provided a very general account of hot 
QCD thermodynamic geometries and systematically ex- 
plored the underlying thermodynamic space as an intrin- 



sic Riemannian manifold, and in particular, introduced 
a covariant geometric notion for both the case of 2-flavor 
and 3-flavor QCD in the realism of quasi-particle theo- 
ries with two distinct chemical potentials. We find in 
both cases that the domain of the thermodynamical in- 
teraction gets shifted towards the origin under the local 
thermal contributions. Hence our geometric construction 
as a covariant intrinsic Riemannian manifold provides a 
realization of the limiting equilibrium thermodynamics 
of such hot QCD near Tc- Our investigation is thus 
an interesting geometric exercise of quasi-particle mod- 
els of QCD, which in turn are described by the quark 
chemical potentials with a transition temperature and 
thus exhibit convincing macroscopic versus microscopic 
duality relations via the fermion number density, as the 
mixing between different flavors/colors encoded in the 
underlying partition function of the corresponding effec- 
tive gauge field theory. In summary, we have presented 
a covariant intrinsic Riemannian geometric perspective 
of any given analytically calculable quark number sus- 
ceptibility tensor, which appears naturally within an ap- 
proximately self-consistent resummation of perturbative 
QCD. It is worth to mention that our investigations thus 
demonstrated are however perturbative in nature, but 
the same line of thought may further be applied for the 
geometric realization of the underlying quark suscepti- 
bility tensors either in the fabric of lattice QCD or that 
of non-perturbative QCD, see for example [2fl. In par- 
ticular, it will be interesting to see whether this kind 
of approach can be pushed further, for example to pre- 
dict geometric properties of chemical correlations, under 
sizable higher order perturbative QCD contributions, as 
well as in general to ascertain geometric features to ex- 
plore the complete non-perturbative QCD. 
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